Exploration of Streptococcus core genome to reveal druggable targets and novel therapeutics against S. pneumoniae.

Streptococcus pneumoniae (S. pneumoniae), the major etiological agent of community-acquired pneumonia (CAP) contributes significantly to the global burden of infectious diseases which is getting resistant day by day. Nearly 30% of the S. pneumoniae genomes encode hypothetical proteins (HPs), and better understandings of these HPs in virulence and pathogenicity plausibly decipher new treatments. Some of the HPs are present across many Streptococcus species, systematic assessment of these unexplored HPs will disclose prospective drug targets. In this study, through a stringent bioinformatics analysis of the core genome and proteome of S. pneumoniae PCS8235, we identified and analyzed 28 HPs that are common in many Streptococcus species and might have a potential role in the virulence or pathogenesis of the bacteria. Functional annotations of the proteins were conducted based on the physicochemical properties, subcellular localization, virulence prediction, protein-protein interactions, and identification of essential genes, to find potentially druggable proteins among 28 HPs. The majority of the HPs are involved in bacterial transcription and translation. Besides, some of them were homologs of enzymes, binding proteins, transporters, and regulators. Protein-protein interactions revealed HP PCS8235_RS05845 made the highest interactions with other HPs and also has TRP structural motif along with virulent and pathogenic properties indicating it has critical cellular functions and might go under unconventional protein secretions. The second highest interacting protein HP PCS8235_RS02595 interacts with the Regulator of chromosomal segregation (RocS) which participates in chromosome segregation and nucleoid protection in S. pneumoniae. In this interacting network, 54% of protein members have virulent properties and 40% contain pathogenic properties. Among them, most of these proteins circulate in the cytoplasmic area and have hydrophilic properties. Finally, molecular docking and dynamics simulation demonstrated that the antimalarial drug Artenimol can act as a drug repurposing candidate against HP PCS8235_RS 04650 of S. pneumoniae. Hence, the present study could aid in drugs against S. pneumoniae.

Microsecond molecular dynamics simulations revealed the inhibitory potency of amiloride analogs against SARS-CoV-2 E viroporin.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) encodes small envelope protein (E) that plays a major role in viral assembly, release, pathogenesis, and host inflammation. Previous studies demonstrated that pyrazine ring containing amiloride analogs inhibit this protein in different types of coronavirus including SARS-CoV-1 small envelope protein E (SARS-CoV-1 E). SARS-CoV-1 E has 93.42% sequence identity with SARS-CoV-2 E and shared a conserved domain NS3/small envelope protein (NS3_envE). Amiloride analog hexamethylene amiloride (HMA) can inhibit SARS-CoV-1 E. Therefore, we performed molecular docking and dynamics simulations to explore whether amiloride analogs are effective in inhibiting SARS-CoV-2 E. To do so, SARS-CoV-1 E and SARS-CoV-2 E proteins were taken as receptors while HMA and 3-amino-5-(azepan-1-yl)-N-(diaminomethylidene)-6-pyrimidin-5-ylpyrazine-2-carboxamide (3A5NP2C) were selected as ligands. Molecular docking simulation showed higher binding affinity scores of HMA and 3A5NP2C for SARS-CoV-2 E than SARS-CoV-1 E. Moreover, HMA and 3A5NP2C engaged more amino acids in SARS-CoV-2 E. Molecular dynamics simulation for 1 µs (1,000 ns) revealed that these ligands could alter the native structure of the proteins and their flexibility. Our study suggests that suitable amiloride analogs might yield a prospective drug against coronavirus disease 2019.

Whole genome sequencing for revealing the point mutations of SARS-CoV-2 genome in Bangladeshi isolates and their structural effects on viral proteins.

Coronavirus disease-19 (COVID-19) caused by SARS-CoV-2 has already killed more than one million people worldwide. Since novel coronavirus is a new virus, mining its genome sequence is of crucial importance for drug/vaccine(s) development. Whole genome sequencing is a helpful tool in identifying genetic changes that occur in a virus when it spreads through the population. In this study, we performed complete genome sequencing of SARS-CoV-2 to unveil the genomic variation and indel, if present. We discovered thirteen (13) mutations in Orf1ab, S and N gene where seven (7) of them turned out to be novel mutations from our sequenced isolate. Besides, we found one (1) insertion and seven (7) deletions from the indel analysis among the 323 Bangladeshi isolates. However, the indel did not show any effect on proteins. Our energy minimization analysis showed both stabilizing and destabilizing impact on viral proteins depending on the mutation. Interestingly, all the variants were located in the binding site of the proteins. Furthermore, drug binding analysis revealed marked difference in interacting residues in mutants when compared to the wild type. Our analysis also suggested that eleven (11) mutations could exert damaging effects on their corresponding protein structures.

Wave-wise comparative genomic study for revealing the complete scenario and dynamic nature of COVID-19 pandemic in Bangladesh.

As the COVID-19 pandemic continues to ravage across the globe and take millions of lives and like many parts of the world, the second wave of the pandemic hit Bangladesh, this study aimed at understanding its causative agent, SARS-CoV-2 at the genomic and proteomic level and provide precious insights about the pathogenesis, evolution, strengths and weaknesses of the virus. As of Mid-June 2021, over 1500 SARS-CoV-2 genomesequences have been deposited in the GISAID database from Bangladesh which were extracted and categorized into two waves. By analyzing these genome sequences, it was discovered that the wave-2 samples had a significantly greater average rate of mutation/sample (30.79%) than the wave-1 samples (12.32%). Wave-2 samples also had a higher frequency of deletion, and transversion events. During the first wave, the GR clade was the most predominant but it was replaced by the GH clade in the latter wave. The B.1.1.25 variant showed the highest frequency in wave-1 while in case of wave-2, the B.1.351.3 variant, was the most common one. A notable presence of the delta variant, which is currently at the center of concern, was also observed. Comparison of the Spike protein found in the reference and the 3 most common lineages found in Bangladesh namely, B.1.1.7, B.1.351, B.1.617 in terms of their ability to form stable complexes with ACE2 receptor revealed that B.1.617 had the potential to be more transmissible than others. Importantly, no indigenous variants have been detected so far which implies that the successful prevention of import of foreign variants can diminish the outbreak in the country.

Genomic Surveillance of SARS-CoV-2 Viruses Collected during the Ending Phase of the First Wave of the COVID-19 Pandemic in Bangladesh.

Mutations, deletions, and the emergence of new variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) may pose a serious health threat. Here, we report the genome sequences of SARS-CoV-2 viruses that were collected from SARS-CoV-2-infected patients during the end phase of the first wave of the COVID-19 pandemic in Dhaka, Bangladesh.

Novel mutations in NSP-1 and PLPro of SARS-CoV-2 NIB-1 genome mount for effective therapeutics.

BACKGROUND: Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), the etiologic agent of coronavirus disease 2019 (COVID-19), is rapidly acquiring new mutations. Analysis of these mutations is necessary for gaining knowledge regarding different aspects of therapeutic development. Previously, we have reported a Sanger method-based genome sequence of a viral isolate named SARS-CoV-2 NIB-1, circulating in Bangladesh. The genome has four novel non-synonymous mutations in V121D, V843F, A889V, and G1691C positions. RESULTS: Using different computational tools, we have found V121D substitution has the potential to destabilize the non-structural protein-1 (NSP-1). NSP-1 inactivates the type-1 interferon-induced antiviral system. Hence, this mutant could be a basis of attenuated vaccines against SARS-CoV-2. V843F, A889V, and G1691C are all located in nonstructural protein-3 (NSP-3). G1691C can decrease the flexibility of the protein. V843F and A889V might change the binding pattern and efficacy of SARS-CoV-2 papain-like protease (PLPro) inhibitor GRL0617. V843F substitution in PLPro was the most prevalent mutation in the clinical samples. This mutation showed a reduced affinity for interferon-stimulated gene-15 protein (ISG-15) and might have an impact on innate immunity and viral spread. However, V843F+A889V double mutant exhibited the same binding affinity as wild type PLPro. A possible reason behind this phenomenon can be that V843F is a conserved residue of PLPro which damaged the protease structure, but A889V, a less conserved residue, presumably neutralized that damage. CONCLUSIONS: Mutants of NSP-1 could provide attenuated vaccines against coronavirus. Also, these mutations of PLPro might be targeted to develop better anti-SARS therapeutics. We hope our study will help to get better insides during the development of attenuated vaccine and PLPro inhibitors.

Recognition of plausible therapeutic agents to combat COVID-19: An omics data based combined approach.

Coronavirus disease-2019 (COVID-19), caused by Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2), has become an immense threat to global public health. In this study, we performed complete genome sequencing of a SARS-CoV-2 isolate. More than 67,000 genome sequences were further inspected from Global Initiative on Sharing All Influenza Data (GISAID). Using several in silico techniques, we proposed prospective therapeutics against this virus. Through meticulous analysis, several conserved and therapeutically suitable regions of SARS-CoV-2 such as RNA-dependent RNA polymerase (RdRp), Spike (S) and Membrane glycoprotein (M) coding genes were selected. Both S and M were chosen for the development of a chimeric vaccine that can generate memory B and T cells. siRNAs were also designed for S and M gene silencing. Moreover, six new drug candidates were suggested that might inhibit the activity of RdRp. Since SARS-CoV-2 and SARS-CoV-1 have 82.30% sequence identity, a Gene Expression Omnibus (GEO) dataset of Severe Acute Respiratory Syndrome (SARS) patients were analyzed. In this analysis, 13 immunoregulatory genes were found that can be used to develop type 1 interferon (IFN) based therapy. The proposed vaccine, siRNAs, drugs and IFN based analysis of this study will accelerate the development of new treatments.

Coding-Complete Genome Sequence of SARS-CoV-2 Isolate from Bangladesh by Sanger Sequencing.

A coding-complete genome sequence of a severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) isolate was revealed. The sample for the virus was isolated from a female patient from Dhaka, Bangladesh, suffering from coronavirus disease-2019 (COVID-19).